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We demonstrate a systematic method for solving the Hamilton-Jacobi equation for general rela- 
tivity with the inclusion of matter fields. The generating functional is expanded in a series of spatial 
gradients. Each term is manifestly invariant under reparameterizations of the spatial coordinates 
("gauge-invariant"). At each order we solve the Hamiltonian constraint using a conformal trans- 
formation of the 3-metric as well as a line integral in superspace. This gives a recursion relation 
~pf+.~, f° r the generating functional which then may be solved to arbitrary order simply by functionally 

■ differentiating previous orders. At fourth order in spatial gradients, we demonstrate solutions for 

' irrotational dust as well as for a scalar field. We explicitly evolve the 3-metric to the same order. 

This method can be used to derive the Zel'dovich approximation for general relativity. 
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In applications to cosmology, Hamilton-Jacobi (HJ) theory plays a vital role in describing the gravitational field. 
It serves as a starting point for a semi-classical analysis used both in stochastic inflation jlj] , Q , || , and in quantum 
C\| , cosmology j|. Recently it has been shown how to apply HJ techniques in solving the classical Einstein's equations. 
For example, one may derive the Zel'dovich approximation describing the formation of pancake structures in a matter- 
dominated Universe |§| , ||] ■ 

The Hamilton-Jacobi equation for pure gravity was first derived by Peres fl. He showed that the momentum 
constraint implies the generating functional or the phase of the semi-classical wavefunctional remains invariant under 
arbitrary reparameterizations of the spatial coordinates (diffcomorphism invariance). However he was not able to 
solve the HJ equation. We find that a crucial ingredient is the introduction of matter fields. We also capitalize on 
the diffeomorphism invariance of the generating functional. 

In this paper we expand the generating functional in a series of spatial gradients. This approximation scheme has a 
long history dating back to Lifshitz and Khalatnikov [Q, as well as Eardley, Liang and Sachs jfj. In a quantum context, 
Teitelboim and Pilati [jllj considered an analogous approach which they called the strong-coupling expansion. 
Accurate to first order in spatial gradients, Salopek and Bond (^] explicitly solved the Hamilton-Jacobi equation and 
the momentum constraint with scalar fields present. In an important generalization, Salopek and Stewart p5§ added 
dust fields and then demonstrated solutions to second order. The principal goal of this paper is to give solutions to 
higher order. 

New techniques yield a practical and systematic method of computing higher order terms. It proves convenient 
to employ a conformal transformation of the 3-metric to simplify the Hamiltonian constraint, which can then be 
integrated using a line integral in superspace. A comparison with an exact solution shows that in situations of interest 
for cosmology even the first few terms can be quite accurate. 

We start in Sec. II by writing out the Hamilton-Jacobi equation for general relativity including two types of matter 
fields, a regular scalar field <j), and a dust field x describing a pressure-free fluid of irrotational massive particles, 
for example, cold-dark-matter. The spatial gradient expansion is introduced, and the Hamilton-Jacobi equation is 
expanded to all orders in spatial gradients. Solutions of this equation correct to zero and two spatial gradients are 
reviewed. A new technique involving a conformal transformation is introduced in Sec. Ill in order to find solutions 
to higher orders. A recursion relation is derived which enables one to compute solutions to any given order, and we 
explicitly give the fourth order terms. For a single dust field, we solve the evolution equation for the 3-metric in Sec. 
IV, and show that our solution recovers the Szekeres exact solution |l3| to Einstein's equations. This is followed by 
some concluding remarks. (Units are chosen so that c = 8ttG = % = 1. The sign conventions of Misner, Thorne and 
Wheeler [Q will be adopted throughout.) 
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II. THE HAMILTON- J ACOBI EQUATION FOR GENERAL RELATIVITY 



In HJ theory the primary object of interest is the generating functional S, which is the phase of the semi-classical 
wavefunctional \P ~ e lS . For general relativity, the HJ equation is simply the energy constraint with the momenta 
replaced by functional derivatives of the 3-metric. In addition we must consider the momentum constraint, which 
demands the gauge invariance of S. 

In this section we write S as a series of terms grouped according to the number of spatial derivatives that they 
contain. The Hamilton- Jacobi equation is likewise written as a sum of such terms. We review solutions for S correct 
to zeroth order (no spatial gradients) and to second order (two spatial gradients) for a single dust field, and for a 
single scalar field. 

The action for general relativity interacting with a scalar field <j> and a dust field x is 



(2.1) 



where ^ R is the Ricci scalar of the space-time metric g^ Vl <j> is a scalar field and x is a velocity potential for irrotational 
dust particles of rest mass m (which is a universal constant). V(<fi,x) is a potential for the matter fields. The rest 
number density n = n(t,x) is a Lagrange multiplier which ensures that the 4-velocity 

= -g^X.v/m (2.2) 

satisfies U^U^ = — 1. x may be rescaled to include to, and so from now on we set to = 1. 
In the ADM formalism the line element is written 

ds 2 = (-N 2 + 1 i3 N l N. j ) dt 2 + 2Nidt dx l + jijdx* dx ] , (2.3) 

where N and JVj are the lapse and shift functions respectively, and 7y is the 3-metric. We can then rewrite the action 
in Hamiltonian form, 

X = J d 4 x (7r*"$ a + ir ij % -NH- N*7U) . (2.4) 

Here the denote the momenta conjugate to 7 y and Q a = (0, x) is an assembly of the scalar and dust fields with 
corresponding conjugate momenta 7r* a . We use the summation convention over the subscript a = 1,2. Variations of 
the lapse and shift functions produce the Hamiltonian and momentum constraints: 

H = 7 ~ 1/ V^ fci [2 7il7ifc - TiiTfei] + y/l + -y ij X,iX,i * x + ^7~ 1/2 {^f + l 1/2 V{<$> a ) 

~\l 1,2 R + ~7 l7 V<Mj = , (2.5a) 
Hi = -2 {fi k ir kj ) tj + n lk 7lk,i + = , (2.5b) 

where R is the Ricci scalar of the 3-metric. 



Variation of the action (2.4) with respect to the momenta give the evolution equations for the field variables 

U - iVV, 4 ) /N = 7 _1 /V, (2.6a) 



(X ~ NXi) /N = yJl + XsX', (2.6b) 
( 7ii - N i{j - N^) /N = 2 7 - 1 /\ kl (2 ljk7il - ll3lkl ) , (2.6c) 



where | denotes covariant differentiation with respect to the 3-metric. Variations of the action ( |2.4| ) with respect to 
the field variables yield the evolution equations for the momenta. However, these remaining evolution equations are 
automatically satisfied provided that 
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satisfy the constr aint e quations (2.5) and provided the evolution equations ( |2.6| ) hold. In t his i nterpretation the 
energy constraint ( 2.5a ) is a self-contained equation for 5, whereas the momentum constraint ( 2.5b ) dema nds that S 
be gauge-invariant. The lapse and shift functions do not enter the constraints. Loosely speaking, equations ( |2.6| ) imply 
that trajectories in superspace are orthogonal to the phase functional. (In a more general formulation, S actually 
generates a canonical transformation to new variables |12|.) 

In terms of the generating functional, the Hamiltonian and momentum constraints are 



H(x) =7~ 1/2 
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Equation (2.8a) is now called the Hamilton- Jacobi equation. 



(2.8a) 
(2.8b) 



A. Spatial Gradient Expansion of the Generating Functional 



Our method will be to expand the generating functional 

5 = 5(0) +5 (2) +5 (4) 



(2.9) 



in a series of terms according to the number of spatial gradients they contain. As a result the Hamilton- Jacobi 
equation can be grouped into terms with an even number of spatial derivatives: 



H = H {Q) +H {2) +W (4) 



(2.10) 



1. Solution of zeroth order Hamiltonian 
The Hamilton- Jacobi equation of order zero is 

= ~S~lij S^ (2 W ,.-7«7«) + ^^ j 
xc(0) 

+ ^ + 7 1 /V($ fl ) = o. 

dx 

The diffeomorphism invariance of the generating functional suggests a solution of the form 

S<°>=-2 J Sx^H^a). 



(2.11) 



(2.12) 



Because the integral is over cfixj 1 ^ 2 the functional is invariant under coordinate transformations, and it clearly 
contains no spatial derivatives. The numerical factor —2 is chosen so that H corresponds to the usual Hubble 
parameter in the long-wavelength approximation. The following condition must be satisfied in order that the zeroth 
energy constraint vanish: 



2 fdHY 2dH 1 



(2.13) 



We now consider two solutions of this nonlinear equation. Firstly for a single dust field and vanishing potential it is 
trivial to show that 
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H{x) = , 
3x 



(2.14) 



is a solution. Setting the shift function to zero, a nd c hoosing our time coordinate to be \ (i- e -> the lapse function 
N = 1), solution of the metric evolution equation (2.6c) gives 



(2.15) 



where kij(x), the 'seed' metric, is an arbitrary function of the spatial variables. This is the standard result for a 
dust-dominated universe at long- wavelengths. 

Secondly we consider a scalar field self-interacting with an exponential potential 



V(<j>) = Foexp 



(2.16) 



where p is a constant. For this case, one can find analytic solutions. The equation (2.13) has an attractor solution 

1 V2 



H(<f>) 
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3(1 - l/3p) 



exp 



(2.17) 



For more details see ref. 



2. Solution of the second order Hamiltonian 

In this section we will give a brief review of the second order, which was treated in Salopek and Stewart [H 
The second order HJ equation is 



H (2) = 2# 7ij - 



8SW dH 5S& 6S& 
2 1 

Sjij d<f> 54> Sx 

1 A<?(°) 1 1 

+p ij X,iX,j~^ ~ ^ R + 2 7 = 0. (2.18) 

Their method was to write an Ansatz for the functional <S' 2 ' which contains all terms that are second order in spatial 
gradients and are diffeomorphism invariant. There are three types of terms — i?, & a i(fr a ' 1 and A$ a . 

5 (2) - J j 1/2 d 3 x (j($ a )R + K afc ($ c )$ a | i $J I + L a ($ b )A$ a ) . (2.19) 

The last term may be integrated by parts, and by redefining 

BT ° 

R ab _^ K ab _ ^_ 20 j 
09b 



we may set L a = 0. Substitution of the Ansatz (2.19) into the Hamilton-Jacobi equation at this order (2.18) yields 
three sets of coupled partial differential equations which then must be solved. The number of equations exceeds 
the number of unknowns. However, no contradiction arises because some are redundant for reasons which will be 
discussed in the next section. 

For a single dust field there are two independent equations. 

QJ 

HJ + — = 1/2, (2.21a) 
ox 



8 J 

2B°4- + K = 
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(2.21b) 



Given the zeroth order solution (2.14), these can be solved to give 
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(2.22) 



The term proportional to the arbitrary integration constant D will be refered to as the complementary functional. 

In an analogous way, for a single scalar field with the zeroth order solution ( 2.17 ), we find (neglecting the comple- 
mentary functional) 



S {2] = J d 3 x7 1/2 



(2.23) 



III. SOLVING THE HIGHER ORDER HAMILTONIAN 

The Ansatz method that was used to solve the second order case can be generalized to higher orders, but it is 
extremely tedious because the number of terms increases dramatically. For example the fourth order case has 15 
types of diffeomorphism invariant terms leading to 15 coupled sets of partial differential equations. Although some of 
the free functions can be set to zero, just as in Sec. II. B. 2, there are still eight variables satisfying fifteen equations. 
It is difficult to proceed further a long these lines, and so we adopt an alternative approach. 



Using the expression for (2.12), we write down the Hamiltonian of order 2n: 



W — - 2^^ + s -^1\ + ti^ - o 



(3.1a) 



Here the remainder term TZ^ 2n \ which is independent of <5>( 2 ™', contains contributions from all previous orders, and is 
assumed to be known: 

„<2„> „,S *S (2p> SS< 2 "-"' „ . ,„^ilSSV'> JS< 2 — 2 r> 

« (2 "> = £ n^-^r (2w " - ™ u) + -< ' g ~, — — 

n /1\ ro(2n-2p) 

+ E(!)(7^X J ) P ^^+V(-), ( 3.l b) 



P =i 

where the superspace potential y( 2 ") is defined to be 



V (2n) = f -i 7 1/2i? + h V V J >,i^ for n = 1, f3 lc) 

[ otherwise. 

Although the H amiltonian constraint of order 2n, (fO), is basically a linear equation for 6>( 2 ™) (unlike the zeroth order 



equation (2.11)), it is nevertheless not im medi ately apparent how to calculate iS^ 2 ™) given the previous orders, partly 



because the term in braces in equation ( 3.1a ) is complicated. In this section we will show, firstly for a single dust 
field, and then for a single scalar field, a technique which yields a recursion relation for the generating functional. 
Finally we indicate briefly how to extend the analysis to multiple fields. 

A. Single Dust Field 



To simplify the HJ equation (3.1) for a single dust field, we utilise a conformal transformation of the 3-metric to 
define a metric variable 

f ij (x)=n- 2 (x(x))- fij (x) (3.2) 

where satisfies 

^=HQ. (3.3) 

dx 



The HJ equation (53a) reduces to 
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(2n) 



(3.4) 



We note immediately that there are some obvious solutions to the homogeneous version of this equation 
SS( 2n ) /Sx\ f =0 which correspond to the complementary functionals. Any functional of only the conformal metric 

1 Jij 

fij that has (2n) spatial gradients and is diffeomorphism invariant will satisfy this equation. For example, at the sec- 



ond order = D J cP xf 1 / 2 R is the complementary functional, where R is the Ricci scalar of /j, 
to the arbitrary part of (2.22). 



This corresponds 



1. The second order revisited 



To show how to use this conformal transformation, we shall now return to the second order HJ equation with one 
dust field present, 



5S^ 2 



in 



= -n^ = \ 1 ^R + 1 ^ d ^x,x, 3 



(3.5) 



The right hand side of this may be expressed in terms of the conformal metric /y using the standard rules for canonical 
transformations, for example, 



r% = n- 2 R\ + n- 1 (n _1 ) . fc f ik - n- 3 f kl n. kl s 



(3.6) 



where R 1 ^ is the Ricci tensor of the 3-metric 7y, R 1 ^ is the Ricci tensor of fij, and a semicolon denotes covariant 
differentiation with respect to /y . Rewriting the right hand side and using the definition of the conformal factor . 
we find 



^ = -K^ [ x (x), f^x)] = /Va (hlR - fi x . Xd ^ H n) - 2/'-' \ : ,..: //<>:') 



5x{x 

By inspection the solution is 

= J d'xf 1 / 2 [j(x)R + k( X )x i X 
where the x-dependent coefficients j and k are 

Ox) 



4fi d X ' + D, k( X ) = Hil , 



(3.7) 



(3.8a) 



(3.8b) 



and D is a constant. We can use the conformal trans formation rules to find the generating functional in terms of the 
3-metric 7^ by using the zeroth order solution (2.14) and the conformal factor 



n = x 2/3 , 



(3.9) 



which satisfies (jyj). This procedure recovers the generating functional ( 2.22 ) obtained by the direct approach. The 
constant of integration D in the definition of j yields the complementary functional at this order. 

However, eq.(3_/7) may be solved directly. A similar situation occurs in potential theory where one wishes to integrate 
an equation of the form 



which has the line-integral solution 



wr gJ{Vk) 



<t>{yk) = J Yl d y'j 9 J (y'i) 



(3.10) 



(3.11) 



G 



where C is an arbitrary contour whose upper end point is yk and whose lower endpoint is arbitrary. If the two 
endpoints are fixed, the solution is independent of the choice of contour C provided that ^2g J dyj is exact. By 
analogy, if / d 3 x d\{x) IZ^ is exact, the solution of eq.(3.7) is 



S (2) 



d 3 x d X '(x) bd(x),fij(x)] 



(3.12) 



The upper limit of integration is x( x ) an( i the lower limit Xo(%) is an independent and arbitrary function. The 
complementary functional is set to zero by choosing xo — 0. We choose the simplest path for the line integral, that 
of a straight line in superspace: 



X'{x)=sx(x), dx'(x) = x(x)ds, < s < 1. 



(3.13) 



The real parameter s is analogous to the Tomonaga-Schwinger proper-time parameter |15| . The integration over s, 
ds s 2 / 3 = 3/5 is trivial leading to 



5 (2) = - Jd 3 x | x ^ (2) bcJij] = J d'xf 1 / 2 ~ X {lnR-fi x ,iX,j§^m) - 2F X &(Hfl) 



(3.14) 



which, after an integration by parts, gives the same result as before, eq.(3.8) 



2. The fourth order 

We illustrate the conformal transformation method for the fourth order. We rewrite IZ^ in terms of the conformal 
metric and we find that, as before, we can integrate the equation directly: it is exact. The resulting generating 
functional written in terms of /y is 



= J d 3 x/ V2 Ux) (& J Rij - j^ 2 ) +m(x)R 2 +n( X ) (& - j X;iX;j + r(x)X-/x ;j X;s 



where 



Kx) 



2 ( X dx> ^ 



Q3( X ') 



m (x) = M , n( X ) 



Kx) 



(3.15a) 



(3.15b) 



The two arbitrary constants, L and M, give the additional complementary functional. Moreover, the complementary 
functional from the second order gives a contribution to this order through the dependence on j. This is a very 
succinct and tidy expression, which becomes more complicated when expressed in terms of the original 3-metric jij, 
(using the conformal factor (|3.9[)) 



^~ 25^^ 3 X\iX\j i+ ggX X\iX y 
where we have neglected all the complementary functionals 



- x *(R«R ij --R^-- X Rx li X 
19 



X X X|iX|j 



(3.16) 



3. The recursion relation 



If for simplicity we choose to ignore all the complementary functionals at each order, we find that a simple recursion 
relation may be derived which relates each order to the previous orders, assuming that the remainder at each order 
is ex act. It is trivial but t edious to show that the remainders 7Z^ and 7Z^ are exact by functionally differentiating 
,(3.8), and , ( |3.15 ). A general proof of the exactness of the remainder TZ^ 2n ^ and, as a result, the arbitrariness 
of the contour C in the integral 



S 



(2n) 



d 3 xd X '(x)K^[x'(x),h(x)} . 



(3.17) 



ic 
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would rely on the fact that the Poisson bracket {H(x), H(y)} yields the momentum constraint: integrability is assured 
provided gauge- invariance is maintained at lower orders. (See also DeWitt and Monc rief w ith Teitelboim jL7|.) 

By induction, one can show that each term in VS 1 ""^ on substitution of the line integral (3.13) has an s-dependence 
of s 2 ™/ 3 . Hence the integral of the HJ equation ( |3.17| ) is 

S (2n) =-J <?xj\s X (x) n^[s X {x)J t] {x)] = J d\ X (x) U^[x(x), f«(x)] , (3.18) 

and we may rewrite this in terms of the 3-metric 7^ 

s (2n) = / d 3 x X (x)n^[x(x)^ 3 {x)\ 



2n + 3 

3 f \ n ~ L 5SW &S( 2 "- 2fe ) 

J d " * ^ g * § (2 W - lijlhl ) 



2n + 3 



+ h ij X,iX tj ) Sx +XV< 2 ») \ . (3.19) 



k=l 



This expression has no reference to the conformal metric /y, and hence is easy to use in practice. In summary, to 
find the generating functional at each order, one must functionally differentiate the previous orders, and substitute in 
the above recursion relation. 



B. One Scalar Field 



To simplify the HJ equation (3.1) for a single scalar field we use a similar conformal transformation to the dust field 
case plus a change of variables for the scalar field. Many of the equations are similar. Once again, fij(x) is defined by 



7y(x) = Vt 2 (u)f i:j (x) 



(3.20) 



where the new variable u(4>) (which at the zeroth order is a comoving, synchronous time variable) is defined through 
an integral 

(3.21) 



)8H 



The conformal factor satisfies the ordinary differential equation 



di 



For a single scalar field, the resulting HJ equation ( 3.1a ) has a form similar to that of dust: 



Su 



n (2n) _ q j 



(3.22) 



(3.23) 



fa 



which again reduces the system to a simple first order functional differential equation with one independent variable. 

W? n > must be rewritten in terms of /y and u. At the second order we find that it is exact, and the HJ equation 
( |3.23 ) has the solution 



S<-V = / d 3 x f 1/2 \j(u)R + k{u)w l u 



where 



j(u)= I -^-Ldu' + D, k{u) = HVL 



(3.24a) 



(3.24b) 



This expression is remarkably similar to th e dust field case ([T^). For the exponential potential considered previously 
(2.16) with the zeroth order solution (2.17), u and the conformal factor are defined by 



8 



u = p/H, Q = u p . (3.25) 



These results can be used to show that the generating functional above ( |3.24 ) agrees with our original approach ( [2.23 ) 



The fourth order solution also has a very similar form to the dust case. 



= J d 3 x f 1/2 J l(u) (sPRy - ^R 2 ^j + m(u)R 2 + n(u) - -/« j U;jttij + r^u^u-j 1 (3.26a) 



where 



l{u) = 2 1 du '9d + L ' m{u) = - / du ' W — + M ' n{u) = ' r{u) = ~4n • 

(3.26b) 



The main difference between this solution and the dust case (3.15) is that here the coefficient of the R term depends 



on u. Expressed in terms of the 3-metric 7^ and <fi this becomes quite complicated. 



2 " 1 (.,,,.. ( 3 1 



1 + Wp^ + (3 - 27) 

The rec ursio n relation for a scalar field may be calculated in a similar way to the dust case. The integral of the HJ 
equation ( [3.2S ) may be written down 



s {2n) = _ d 3 x du f^ TZ ( - 2n) [u'(x),f ij {x)} , (3.28) 



where uo is a constant function which we can set to zero (this sets the complementary functionals to zero). Once 
again, as in the dust case (Sec. III. A. 3), the exactness of the second and fourth order remainders can be shown 
explicitly. The exactness of the higher orders would seem to follow from the Poisson brackets of the constraints. One 



can solve ( p. 28] ) using a line integral as before 

u'(x)=su(x), du'{x) =u(x)ds, < s < 1. (3.29) 

By induction, one can show that each term of US 2 ™} , on substitution of this line integral, has an s-dependence 

s (2n-2) + (3-2«)p and hence 



n-1 



1 /• 3 p [ ^ ^ SS^US^ 2n -^ 



d3x m^ 7 " 1/2 E JZ— - 7«7«) 



(2n - 1) + (3 - 2n)p J H(4>) [ 1 £j £ 7ij 5 lM 

Again, this expression has no explicit reference to the conformal metric, and is easy to use. 

C. Multiple Fields 

One can solve the Hamilton- Jacobi equations for multiple fields by employing the method of characteristics. The 
characteristic equations are 
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d\ d<f> 



with solution 



(3.31a) 



Itj = ^ 2 (x(x)) fij(x) 



(3.32a) 



(f) = h(x(x) + q(x)), 

where the fij(x) and the q(x) are constants of integration along a characteristic, and 



an 

dx 



HCl. 



In this way equation (3.1a) has only one independent field X (x 



D X (x) 



n (2n) 



fij ,1 



(3.32b) 



(3.33) 



(3.34) 



where D denotes functional differentiation along a characteristic. One of the characteristic equations (3.32a) can be 
seen to define a conformal transformation. 



IV. THE EVOLUTION OF THE 3-METRIC 



Using an iteration scheme, we show how to solve the evolution equation (2.6c) for the 3-metric. We give explicit 
solutions to fifth order for a homogeneous dust field. A comparison is made with the Szekeres exact solution of 
Einstein's equations. The ZePdovich approximation is briefly discussed. Firstly we find it useful to define a restricted 
generating functional. 



A. The restricted generating functional 



At the end of the calculation, we are free to choose our time hypersurface. For a single dust field, a natural choic e 
is a surface of uniform Xi which corresponds to comoving synchronous gauge — the evolution equation for x-> ( 2.6f| ), 
shows that the lapse N = 1. For now on, we will also choose the shift function N l = 0. 

We define the restricted generating functional, S [jij{x)\x], in which the dust field is assumed to be homogeneous, 
by 



^ b/ij(x)\ X ] = S [jij(x), x(x) = X ] ■ 
The rate of change of the restricted generating functional <S[7ij(x)|x] is 



dS[^j(x)\x} 
dx 



= d 



SS['yi j (x) ) x(x)] 



Sx(x) 



xO)=x 



in which case the Hamilton- Jacobi equation implies the integrated Hamilton- Jacobi equation for iSfyy (x)|x] 



as 

dx 



+ <fx{-i 



-1/2 



and the momentum constraint reduces to 



= -2 T: 



5S 



),k- 



ss 



$lkj(x) J j S^i k (x) 



(4.1) 



(4.2) 



(4.3) 



(4.4) 
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The restricted generating functio nal i s much easier to calculate in practice. The gradient expansion and the conformal 
transformation trick work for eq^O) as well. However, integrability questions are not important here because there 
is only one equation to solve which always possesses solutions. 



We can either compute the restricted functional using the integrated HJ equation (4.3), or using the full functional 
(3.16) and setting x,i = at the end of the calc ulati on. Both methods give the same answer. The restricted generating 
functional to fourth order is easily seen from ( 3.16 ). 



-.„(.<-)|\] = - /^ 1/2 ^x 3 



Ha — — R 2 



(4.5) 



One can go even further and compute the sixth order functional for a single homogeneous dust field: 

S {6) hijWlx] = J d 3 x^ 2 1 {-^i? 3 + Y RRijR * -±R lJ R\Rjk + R lJ R m ]k - \RR ]k \u} ■ (4.6) 



B. The evolution of the 3- metric for dust 

One can use these results to compute the evolution equation for the three-metric correct to higher order in spatial 
gradients for a single dust field. We a dopt the iteration scheme proposed in reference |12|. The evolution equation 
computed from the ADM action (2.6c) on substitution of the spatial gradient expansion of the generating functional, 
truncating at sixth order is 



ox V o^fki 



Sjki S"f k i Sj k i J 



(4.7) 



where we have set the shift Ni to be zero. We choose our time coordinate to be \ which fixes the lapse function 
N = 1. 



Retaining only the leading term on the right-hand side gives the first order result for 7^ 1 ^(x,a;), eq.(2.15). Substi- 
tuting this expression into the right-hand side, and retaining only third order terms, one finds 



(3) 



(4.8) 



where Rij = Rij{ki m ) and R = i?(fc; m ) are the Ricci tensor and Ricci scalar, respectively, of the seed metric ki m . The 
third-order expression for the 3-metric is then: 



l\V (X,x) = X 4/3 h 3 (x) + -x 2 R(x) k l3 (x) - AR i3 (x) 



20' 



(4.9) 



One more iteration leads to the fifth order expression 



„(5) 



= x A/3 hj + ^x 2 (Rk ij -iR ij ) 



81 



^x 8/3 ( { ^R z + z R 



20' 



-10RR; 



-R-ij + HR^Rjn - -Rij-,k' k 



350' 
5 

— j 

2 



89 



32 



4R kl R k i } k 



(4.10) 



where a semi-colon denotes covariant differentiation with respect to fcjj . A similar iteration scheme could be applied 
equally well to the scalar field case. 

1. Comparison to an exact solution 

In this section we will compare our approximation scheme to an exact solution of Einstein's equations — the 
axisymmetric specialization of the Szekeres solution [n3[ for irrotational dust, which has a line element 



ds 2 = -dt 2 + a 2 {t) [dx 2 + dy 2 + (d(x l ) - c{z)a(t)) 2 dz 2 



(4.11) 
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where d(x l ) has the form 



d{x l ) = d (z) ~ -c(z)(x 2 +y 2 ), 

and do(z) and c(z) are arbitrary functi ons of z. 

We interpret our evolution equation ( 4.10 ) as evolving the seed metric kijix) which we choose to be 



diag [l,l,d 2 (x i )] 



(4.12) 



(4.13) 



where d(x % ) is given by (4.12). We can now evaluate all the terms on the right hand side of (4.10). It is w orthw hile 
using a computer algebra program, particularly for the fourth order part. Even for a simple metric such as ( 4.13| ) the 
individual terms such as Rij-k' k are quite complicated, and so it is extraordinary that the end result is so simple. We 
find that 



7 | 7 5) = x 4/3 diag [l,l,(d(^)-c(z) x 2 / 3 ) 2 



(4.14) 



which is precisely the Szekeres solution ( 4.11 ) ! 

We can use this to check that there is no seventh order contribution to the Szekeres soluti on fr om our method by 
finding the solution to the sixth order evolution equation. We find that with the seed metric (4.13) there is indeed no 

contribution to the metric from the sixth order, i.e., j^J = 7}? '. 

For situations of interest for cosmology, we do not wish to assume such a special form for the seed metric. For this 
case, the expansion can be improved quite significantly if one notes that the the deter min ant of the 3-metric must 
be non-negative. However, one disturbing property about the third order expression (4.9) is that after a sufficient 

amount of time, the 3-metric may no longer be positive definate not only doe s the expansion break down, it gives 

nonsensical results. One may remove this embarrassing problem by writing (4.9) as a 'square': 



7? = X 4/3 



hi 



9 2/3 

40 X 



Rk a - AIL, 



2/3 



Rk 



4R 



jm 



(4.15) 



Both eqs.(4.S) and ( 4.15 ) agree to third order, but the latter is the preferred form because it is non- negative. This 
improved result reproduces the Zel'dovich approximation. See ref. Q and |l{| for a discussion of higher order terms. 



V. CONCLUSIONS 



In this paper we have developed a systematic method for solving the Hamilton-Jacobi equation for general relativity 
using a spatial gradient expansion. We derived a recursion relation by which we can compute the generating functional 
at each order by summing various combinations of the functional derivatives of lower orders. The key ingredients are 
a conformal transformation of the 3-metric and a line integral in superspace. 

Our formalism provides some deep insights into the structure of semiclassical superspace which now far exceeds 
investigations in homogeneous models. Superspace describes an ensemble of evolving universes, and its complexity 
strains the imagination. However, the gradient expansion allows one to separate superspace into an infinite but 
discrete number of manageable pieces which are relat ively easy to understand. 

Because the line integral in superspace (see eq.( 3.12 )) is independent of the choice of contour (provided the endpoints 
are fixed) , we now begin to understand the invariance of the generating functional under different time-hypersurface 
choices. The exactness of the remainder terms is highly nontrivial, and this important property can essentially be 
traced to the fact that the Poisson bracket of the Hamiltonian constraints, {H(x),'H(y)}, yields the momentum 
constraint. Integrability is closely related to the gauge-invariance of the theory. 

The first few terms of the gradient expansion are of considerable interest for a dust-dominated Universe describing, 
for example, the cold-dark-matter model. It is quite reassuring that the fifth order 3-metric reproduces the exact 
Szekeres solution. In addition, HJ theory provides an efficient and practical means of calculating higher terms in 
the Zel'dovich approximation ||, |19|. However, we have not discussed the convergence of the expansion. For more 
general problems, it seems likely that one must employ techniques that will improve the rate of convergence pM. 
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